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ABSTRACT 

A  stabilization  problem  for  Burgers’  equation  is  considered.  Using  linearization,  vari¬ 
ous  controllers  are  constructed  which  minimize  certain  weighted  energy  functionals.  These 
controllers  produce  the  desired  degree  of  stability  for  the  closed-loop  nonlinear  system.  A 
numerical  scheme  for  computing  the  feedback  gain  functional  is  developed  and  several  nu¬ 
merical  experiments  are  performed  to  show  the  theoretical  results. 
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1.  Introduction 


Burgers’  equation  is  a  one-dimensional  simple  model  for  convection-diffusion  phenomena 
such  as  shock  waves,  supersonic  flow  about  airfoils,  traffic  flows,  acoustic  transmission,  etc. 
For  high  Reynolds  number,  the  solution  of  Burgers’  equation  produces  steep  gradients  due 
to  the  nonlinear  nature  of  the  convection. 

In  this  paper,  we  consider  a  stabilization  problem  for  Burgers’  equation.  We  employ 
“a-shifted”  linear  feedback  control  laws  to  obtain  a  desired  degree  of  stability,  on  a  certain 
energy  space,  for  the  closed-loop  nonlinear  system. 

Burgers’  equation 

d  d  <92 

(i-1)  dtz^'x^ +  = 

was  introduced  by  Burgers  [5,6,7]  as  a  simple  model  for  turbulence,  where  e  >  0  is  a  viscosity 
coefficient.  Since  then,  many  researchers  have  considered  the  conservation  law 

d  d 

(1.2)  ~z(f,x)  -f  z(t,x)— z(t,x)  -  0 
and  the  “viscosity  solution” 

(1.3)  z(t,  x)  =  lim  zl(t,  x), 

<10 

where  z'(t,x)  satisfies  equation  (1.1),  see  [9,16,18,22,24,25,26]. 

Oleinik  [26]  proved  that  for  any  Z/°°-initial  data,  there  is  a  unique  viscosity  solution  for 
equation  (1.2)  and  the  solution  satisfies  the  “entropy  condition” 

z(t,x  +  a)  -  z(t,x)  E 

(1.4)  -  <  — 

a  t 

for  all  t  >  0,  n  >  0,  — oo  <  x  <  oo  and  for  some  constant  E  >  0.  A  complete  discussion  of 
these  results  may  be  found  in  [32]. 

Almost  no  results  exist  for  the  control  problem  associated  with  Burgers’  equation.  Chen, 
Wang  and  Weerakoon  [8,35]  considered  an  optimal  control  problem  for  equation  (1.2)  with 
— oo  <  x  <  oo.  The  problem  was  to  select  an  initial  function  to  minimize  a  specific  cost 
functional  J.  They  obtained  sufficient  conditions  for  the  differentiability  of  J  with  respect  to 
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the  initial  function  and  explicit  expression  of  the  entropy  solution  of  (1.2)  in  terms  of  initial 
data. 

In  this  paper  we  consider  a  control  problem  for  Burgers’  equation  (1.1)  defined  on  a 
finite  interval.  Specifically,  we  will  find  several  feedback  laws  stabilizing  the  nonlinear  system 
(1.1)  with  a  certain  exponential  decay  rate.  The  feedback  laws  will  be  obtained  from  the 
linearized  equation.  Curtain  [1 1]  has  considered  a  stabilization  problem  for  certain  semilinear 
evolution  equations.  Using  Kielhofer’s  stability  results  for  semilinear  evolution  equations  [20], 
she  showed  that,  under  certain  conditions,  there  exists  a  finite  dimensional  compensator 
which  produces  a  stable  closed-loop  system.  These  finite  dimensional  compensators  are  also 
obtained  from  the  linearized  control  system.  Applying  her  results  to  Burgers’  equation  (1.1) 
with,  for  example,  Dirichlet  boundary  conditions,  one  can  obtain  the  stabilizability  results 
of  the  closed-loop  system  which  are  similar  to  ours.  However,  in  [11],  there  is  a  restriction  on 
the  action  of  the  output  operators.  The  domain  of  the  output  operator  was  required  to  be 
a  certain  subspace  of  L 2  which  contains  the  Sobolev  space  H^.  In  this  paper,  we  investigate 
optimal  feedback  laws,  instead  of  finite  dimensional  compensators,  in  the  sense  that  they 
minimize  certain  energy  functionals. 

Well-posedness  and  stability  results  for  the  open-loop  system  are  obtained  in  Section  2. 
In  Section  3,  a  “shifted”  linear  control  problem,  (LQR)a,  is  introduced.  Under  appropriate 
selection  of  input  and  output  operators,  (LQR)„  is  stabilizable  and  detectable.  The  feedback 
control  law  obtained  from  (LQR)q  produces  the  desired  degree  of  stability  for  the  closed- 
loop  nonlinear  system,  (Theorem  3.10).  In  Section  4,  a  numerical  scheme  for  computing  the 
“feedback  functional  gains”  is  developed  and  several  numerical  experiments  are  performed. 

We  shall  use  standard  notation.  If  (X,  jj  •  ]]*)  and  (Y,  [|  •  ||r)  are  normed  linear  spaces, 
then  C(X ,  Y)  will  denote  the  space  of  all  bounded  linear  operators  from  X  to  Y  and  for  any 
A  €  C(X,Y),  ||A||  or  ||/4||c(x,k)  will  denote  the  operator  norm  on  the  space  C{X,  Y).  In 
the  event  that  X  =  Y  we  denote  C{X,Y)  by  C(X).  From  time  to  lime  we  will  use  ||  •  || 
without  any  subindex  for  vector  or  operator  norm.  In  all  such  cases  the  appropriate  index 
for  ||  •  ||  will  be  understood  from  the  context.  For  a  Hilbert  space  X,  we  denote  the  inner 
product  on  X  x  X  by  <  •,•  >x ■  Given  a  linear  operator  A  from  A'  into  itself,  we  denote 
its  domain,  spectrum,  resolvent  and  adjoint  by  V(A ),  cr(A),  p(A)  and  A*,  respectively.  For 
real  numbers  a,  b  with  a  <  b,  Lp(a,b\  X),  1  <  p  <  oo,  will  be  the  space  of  all  Lebesgue 
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measurable  functions  /  from  (a,b)  to  X  such  that  ||/||z,p(a,6)  =  (/a  \f(x)\pdz)?  <  oo.  The 
spaces  Hk(a,b)  and  //o(a,6)  are  the  standard  Sobolev  spaces  defined  by  Hk(a,b)  =  {/  6 
L2(a,6)|/(;>  €  L\a,b),j  =  0,1,  ••  ,*}  and  //*(«, 6)  =  {/  G  7/*(a,  6)|/^(a)  =  /<»(&)  = 
0,j  =  0, 1 ,  -  •  -  ,  k—  1},  respectively.  The  dual  space  H~k(a,b)  of  //g(a,6)  is  the  space  of  all 
continuous  linear  functionals  on  H^a^b)  represented  by  the  inner  product  <  >Li(a,b)- 

Finally,  we  present  a  physical  example  that  may  be  found  in  most  standard  references  to 
motivate  the  control  problem  for  Burgers’  equation.  Other  examples  involving  traffic  flows, 
supersonic  flow  about  airfoils,  acoustic  transmission  and  turbulence  in  hydrodynamic  flows 
can  be  found  in  [12]  and  the  references  given  there.  The  following  example  is  taken  almost 
directly  from  [12]. 

Example  (Shock  Waves) 

An  impulsively-started  piston  moving  at  a  constant  velocity  into  a  tube  containing  a 
compressible  fluid  initially  at  rest  creates  compression  waves.  The  compression  waves  even¬ 
tually  coalesce,  due  to  the  nonlinear  nature  of  the  convection,  to  form  a  single  shock  wave. 
The  one-dimensional  unsteady  motion  of  the  fluid  is  governed  by  the  continuity  equation 

d  d  $ 

(1.5)  p(t,x)  +  p(t,x)~  v{t,x)  +  v(t,x)-^p(t,x)  =  0 

and  the  x-momentum  equation 

c)  d  d  ^2 

(1.6)  dtv(t,x^  + 

where  p  is  the  density,  v  is  the  velocity,  p  is  the  pressure  and  6  is  the  “diffusivity  of  sound”. 
It  is  convenient  to  replace  the  density  by  the  sound  speed,  a  =  a(t,x )  by  a(Ax)/a0  = 
(p{t,  x)/p0)2^~,  where  7  >  1  is  the  specific  heats  ratio  and  the  subscript  0  refers  to  the 
undisturbed  values  [23].  Equations  (1.5)  and  (1.6)  become 

d  d  7  —  1  d 

(1.7)  dt,a('t'X'>  +  v(t'x)g^a(t^x)  +  — y~  a{t,x)—v{t,x)  =  0 
and 

d  d  2  d  d2 

(1.8)  dtv(l'x )  +  v(<,x)^t,(i, J)  + 
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where  8  is  a  function  of  the  undisturbed  (to  the  right  of  the  shock)  viscosity,  density,  specific 
heat  and  thermal  conductivity  of  the  medium.  Equations  (1.7)  and  (1.8)  can  be  simplified 
by  introducing  the  Riemann  invariants, 

a(t,x)  v(t.x)  a(t,x )  v(t,x) 

r(t,x)  =  - r  + —^— ,  -<^x)  = - ; - 5— 

7  —  1  2  7  —  1  2 

to  give 

<7  «7  r 

(1.10)  ^r(l,i)  +  (a(t,x)  +  v(t,x))—r(t,x)  =  ^-^ir{^x)  “  si^x)) 

and 

(Lii)  dts^'x^  ~  ^ t,x ^  ~  v(^x))^s(^x)  =  7^x r(-s(h  j)  -  r(<,i)). 

Consider  the  propagation  of  a  disturbance  into  an  initially  undisturbed  region,  s  =  ,s0 
where  ,s0  =  Then  the  problem  is  governed  by  equation  (1.10).  Hut  from  equation  (1.9), 

(1.12)  a(t,x)  +  v(t,x)  =  ~^-r(t,x)  + 

thus  equation  (1.10)  becomes 

d  7+1  7  —  35  8  d2 

(U3)  -r[t,x)  +  (— r  (i.x)  +  — )^r((..r)  =  jgjjr((,.r). 

As  the  final  step  we  introduce  the  change  of  variables 

(1.14)  z{t,x)  =  1^-{r(t,x)  -  r0),  £  =  x  -  a0t 

to  give  Burgers’  equation 

(1.15)  !*«,<)  +  =(<,oJC(<,{)  =  —*(»,{). 

From  equations  (1.12)  and  (1.14)  we  have  z(t,£)  =  {n(/,£)  +  e(/,£)}  —  {  e0  +  Oy},  wliere 
t’o  =  0,  i.e.,  m(<,0  is  the  excess  wavelet  velocity  (the  difference  between  propagation  speeds 
of  disturbance  in  stagnation  and  nonstagnation  conditions).  The  coordinate  £  is  measured 
relative  to  a  frame  of  reference  moving  with  the  undisturbed  speed  of  sound  n0. 
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2.  Well-posedness  and  stability  of  Burgers’  equation 

In  this  section,  we  consider  well-posedness  and  stability  properties  of  the  solution  for 
Burgers’  equation  with  Uirichlet  boundary  condition.  These  results  will  be  needed  in  the 
analyses  of  our  control  problems  in  the  next  section.  We  first  consider  an  abstract  version 
of  this  problem  and  then  specialize  to  Burgers’  ecpiation. 

Consider  an  initial  value  problem 

(2.1)  —  z(t)  =  Az[t)  +  f(t.,z(t)),  z(l0)  =  z0,  (t  >  t0), 

at 

on  a  Hilbert  space  X ,  where  A  is  the  infinitesimal  generator  of  an  analytic  semigroup  S(t ) 
satisfying  ||5(t)||£(X)  <  Me"'1,  t  >  t0 ,  for  some  constants  M  =  M( ui)  >  1  and  u>  >  0.  Since 
S(t)  is  analytic,  the  fractional  powers  of  A\  =  —  A  +  al  are  well-defined  for  any  a  >  u  [27, 
Chapters  1,2].  Since  0  £  p{Ai),  the  resolvent  of  A\,  A ^  is  invertible  for  all  0  <  p  <  1. 
Therefore,  the  graph  norm  ||z||  -f  ||^z||  on  the  domain  V(Ay)  of  Af  is  equivalent  to  the 
norm  ||z||^  =  ||^4^z||.  VVe  denote  the  Hilbert  space  V(A^)  with  the  norm  ||z||  -f  1 s 1 1  or 
MM|by(X„||-lU. 

We  shall  make  the  following  assumption  [17,27]. 

Hypothesis  (F):  Let  U  be  an  open  subset  of  [?0,oo)  x  .  The  function  /  :  U  — ♦  X 
satisfies  the  hypothesis  (F)  if  for  every  ( t,z )  6  U  there  is  a  neighborhood  V  C  V  and 
constants  L  >  0,  0  0  <  1  such  that 

(2-2)  ||/(Hi  z\)  ~  f{t 2’>zi')\\x  ^  L{\ti  —  t2  f  +  j|zj  —  z2 ||M) 

for  all  ( ,  zr, )  6  V,  i  =  1,2,  i.e.,  /  is  locally  Holder  continuous  in  /,  locally  Lipschitzian  in  z, 
on  U. 

Now  we  are  ready  to  state  the  local  existence  theorem  for  the  solution  of  equation  (2.1). 
The  following  result  appears  as  Theorem  3.3.3  in  [17]. 

Theorem  2.1.  Let  A  be  as  before  and  /  satisfy  hypothesis  (F).  Then  for  any  (to,~o)  € 
U  C  R+  x  Xp,,  there  exists  T  =  T(<0,z0)  >  0  such  that  equation  (2.1)  has  a  unique  (strong) 
solution  z{t)  on  [<0,  +  T)  with  initial  value  z(t0)  =  z0. 

Now,  consider  Burgers’  equation,  with  Dirichiet  boundary  condition,  on  a  finite  interval 
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M  given  by 

r\  r\ 2  r\ 

— z(t,x)  =  e—  z{t,x)  -  z{t,x)—z{t,x),  0  <  x  <  (,  t>  0 

z(t,0)  =  z(tj)  =  0 
z(0,x)  =  z0(x), 

where  e  =  ~  >  0  and  Re  is  the  Reynolds  number.  In  order  to  place  the  system  (2.3)  into 
a  semigroup  framework  let  z(t)(-)  =  2o(‘)  =  2(05-)  and  H  =  L2( 0,£).  Define  the 

operator  At  by 

(2.4)  At<j)  =  t  <t>" 

for  all  <f>  G  T>(AC)  =  H2(0,l)  H  Hq( 0,£).  The  system  (2.3)  can  now  be  written  as  the  initial 
value  problem 

(2.5)  4 2(t)  =  Atz(t)  +  /(*, 2(0)>  «(0)  =  2o,  (t>  0), 

at 

on  the  space  //,  where  f{t,z)  =  —  zz'  is  defined  on  the  space  It  is  well-known 

[17,27,34]  that  A(  generates  an  analytic  semigroup  S(t)  on  H. 

We  summarize  the  basic  properties  of  the  infinitesimal  generator  A(  and  its  semigroup 
S(t),  t  >  0,  in  the  following  remark. 

Remark  2.2.  (i)  The  spectrum  a(A()  of  Ae  consists  of  all  eigenvalues  An  =  -en27r2/£2, 
n  =  1,2,  •  •  • ,  and  for  each  eigenvalue  An  the  corresponding  eigenfunction  </>„  is  given  by 

(2.6)  <f>n(x)  =  V2  sin  —  x,  0  <  x  <  L 

(ii)  The  operator  A(  is  self-adjoint,  i.e.,  A<  =  A%  and  the  semigroup  S(t)  can  be  repre¬ 
sented  by  the  following  formula 

OO 

(2.7)  5(02  =  X]  <z,K>  4>n 

n=  1 

for  all  z  £  H,  where  <f>n' s  are  defined  by  equation  (2.6).  Moreover,  from  equation  (2.7),  it  is 
easy  to  see  that  S(t)  has  the  stability  property 

(2.8)  ||S(I)|U,h,  <  e'"’’"'",  <>0. 
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A  simple  application  of  Schwartz  inequality  gives  the  following  first  Poincare  inequality 
[36,  p.  116]. 

Lemma  2.3.  For  any  2  G  77,5(0, 7), 

(2-9)  IMIh<«Mh, 

where  77  =  75(0,7). 

Remark  2.4.  (i)  The  above  lemma  gives  an  equivalent  norm  \\z\\Hi  =  \\z'\\i2  on  the  space 

#<5(0,7). 

(ii)  It  is  well-known  that  V((—A()%)  =  77,5(0,7)  [17,  p.  29], [21,  p.  326]. 

Lemma  2.5.  For  any  2  G  77,5(0,7)  =  £>((  —  j4£)5),  the  following  inequalities  hold. 

(2.10)  I|5(0*IIh>  <e“‘ylIWIn»,  (<>0) 

(2.11)  l|S((WI„<  (*>0). 

where  7  =  £7r2/72. 

Proof:  For  any  2  G  77,5(0, 7),  we  know  that  (-A£)s5(0*  =  S(t)(-At)* 2.  Hence,  by 
Remark  2.2, 

WMh*  =  \\S(t)z\\L2  +  \\(-A^S(t)z\\L2  <  ||5(0||  ( INI  +  ll(-4£)Ml  )  <  e-7*  INI//- 

The  inequality  (2.11)  follows  from  Remark  2.2,  Lemma  2.3  and  the  estimate 

OO 

I|S(<WIh.  <  (1  +  0  ||(-/l,)iS(l)Z||„  =  (1  +  t)  ||(-/t,)*  Yi  «*"'  <  Z,*n  >  *,|| 

n=  1 
OO 

=  (1  +0  II  ~TeKt  <  Z^n>  V^COS  ^—~x  \  | 

71=1 

<  (1  +  7)  (sup{~e(An+7)‘  :  n  =  1, 2, . . .  })e_7t|N|// 


sup(2Ie',-+’)':n=l,2, ••■}<( 

1  If.  <>&• 


where  An  =  —en27r2/£2  and  <j>n(x)  =  \/2  sin  ^x,  n  =  1 , 2,  •  ■  • . 


Remark  2.6.  The  inequality  (2.11)  holds  for  every  :£//  =  Z/2(0,£),  since  the  semigroup 
S(t)  is  analytic. 

Now  we  have  the  well-posedness  and  stability  properties  of  Burgers’  equation  (2.3)  on 
the  space  Hq( 0,f).  The  following  theorem  follows  from  an  application  of  Theorem  5.1.1  in 

[U]. 


Theorem  2.7.  For  any  given  /3  >  0,  0  <  /?  <  7  =  £tt2/£2,  there  is  a  p  —  p((,e,/3)  >  0 
such  that  for  any  initial  data  Zo  G  Hq(Q,£),  with  1 1  ~  o  1 1  5;  §’  there  is  a  unique  solution 

z(t)  =  2(/0;2O)  G  //o(0,f)  of  equation  (2.3).  Moreover,  the  solution  satisfies  the  inequality 


(2.12) 


r(t,0;  20)||//i  <  2e  ^||z0||//i  (t  >  0) 


and  p  =  p(i,  e,fi)  >  0  can  be  chosen  to  satisfy 


(2.13) 


0  <  p  < 


>/^(7-/*) 


v/27T  (1  +  +  V/27T£  ) 


Proof:  Note  that  0  is  an  equilibrium  point  for  the  system  (2.5).  Since  X>((  —  At)* )  =  Hq, 
if  the  nonlinear  term  f(z)  =  —  zz'  satisfies  the  hypothesis  (F)  with  index  p  =  then,  by 
Theorem  2.1,  we  have  a  unique  local  solution  z(t,0;2O)  on  the  space  //q.  It  is  easy  to  see 
that  ||/(2,)  -  f(z2) \\L7  <  (|ki||;/oi  +  ||z2||//i)lki  -  z*\\hI  for  all  zuz2  G  H£,  uniformly  in 
t  >  0.  Hence.  /  satisfies  the  hypothesis  (F). 

For  the  global  existence  and  uniqueness  of  the  solution  2(t,0;2o)  G  Hq,  let  Zo  be  any 
initial  data  in  //0*  with  I  l^o  |  I//q  <  |,  where  p  =  p((,e,/3)  satisfies  the  condition  (2.13).  It 
follows  that 


(2.14) 


ds  + 


(1  +  Qtt  ^1 

f(7  -  /?)  '  2  ’ 


since  JJ°  e  ^  ds  = 


— )^=,  where  F  is  the  Gamma  function  and  T(  h)  = 

V2/ 


o 1 1 //^  <  Then,  by  the  local  existence  property,  there  is  a  unique  solution 


/-,-/?  rG)sinf 

x/ir.  Let 

z(<,0;zo)  G  //q  satisfying  the  inequality  ||z(f,0;  z0)||/^  <  p  on  an  interval  [0,t,)  for  some 
t\  >  0,  where  <,  is  chosen  as  large  as  possible.  We  will  show  that  t,  must  be  infinity.  Suppose 
that  t\  is  finite.  Then  we  must  have  ||z(t  .)lks  >  p.  Note  that,  on  the  interval  [0,/,), 


(2.15) 


ll/M0)llf  =  II  -  .'(<).-'(<)lly  <  (IWOIM2  <  p '. 
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whore  '  -  4~.  Lemma  2.5,  Remark  2.6  together  with  inequalities  (2.14),  (2.15)  yield 


IMMIk  =  ||S(M*0+  P  Ik 

Jo 

<  l|S(M*bllw«  +  P  11^1  -  s)M*)) ll«*  * 

Jo 


-5*-i 


This  is  a  contradiction.  Therefore,  the  unique  global  solution  c(L0;-o)  exists.  Moreover, 
from  the  above  estimate,  we  know  that  if  H^e';//'  5:  §  then  ||r(L  0;  ~o)|]//i  <  P  f°r 
t  €  [0,  oo). 

Finally,  we  will  derive  the  stability  result  (2.12).  Let  w(t)  =  sup{  1 1 ■? ( ^ ) 1 1 // » t03  ■  0  < 
s  <  l  }.  We  then  have 


IWOlUic'”  <<•'"(  l|S(0*-oll„,  + 


fm  -»)/(*-(» 

Jo 


I//' 


cLs) 


<  Ikollm  + 


i+f 


SriJ 


VT('  Vo  \/T 


:e-(y~^(t-)ds)w(/)  + 


(l+f);r 


’-/,(/  c-^-^ds  )w(t) 

Jo 


<  ll2o||//, 


■.  +/>{  -Ltf  r -Lc-o-i»- ,l,  +  <‘  +  )  „,(D  <  ||.-0||„,  +  !.,.(()• 

0  io  W  ^(7  -  P)  0  2 


Therefore,  ?e(<)  <  2 1 1  ^r0 1 1  /f  *  and  ll-(0li/^  <  2c  ^H-oll/^- 

Remark  2.8.  Rankin  [29]  considered  well-poscdness  properties  for  a  certain  type  of  semi- 
linear  evolution  equations  where  the  nonlinear  terms  are  in  divergence  form.  According  to 
his  results,  we  can  see  that  equation  (2.3)  has  a  unique  (strong)  solution  for  initial  data  in 
fp(Qj),  p  >  4.  To  get  this  result,  he  used  the  analyticity  of  the  semigroup  S(t)  and  the  fact 
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that  the  differential  operator  —  on  Hq(0,  I)  can  be  represented  by 


(2.16)  =  B 

for  some  bounded  operator  B  E  £(//),  where  II  =  L2(0,I).  In  general,  it  is  not  true  that 
s/e  ~  =  (  —  A()5 .  This  result  could  be  used  to  analyze  the  stability  property  of  the  solution 
of  Burgers’  equation  with  initial  data  in  Lp(0,().  p  >  4. 

3.  Linear  control  problem 

As  we  noted  in  Section  2,  the  open-loop  (no  control)  solution  of  Burgers’  equation  (2.3) 
decays  exponentially  in  the  topology  of  the  energy  space  H o(0,f)  (see  Theorem  2.7).  How¬ 
ever,  the  decay  rate  depends  on  the  viscosity  e  >  0.  We  now  explore  the  possibility  of 
obtaining  an  exponential  decay  rate  independent  of  viscosity  by  feedback  laws. 

The  basic  model  is  governed  by  an  abstract  system  of  the  form 

^  ^  ^z(t)  =  Az{t)  +  Bu(t),  z(0)  =  zoeH, 

y(t)  —  Cz(t),  t  >  o, 

where  H,  U  and  Y  are  Hilbert  spaces,  u(-)  £  L2{0,oc;U),  y(-)  E  Z/2(0,oo;y),  and  A  is  the 
infinitesimal  generator  of  an  analytic  semigroup  S(t)  on  a  Hilbert  space  II.  Assume  that 
B  E  C(U,  //),  C  E  C{H,Y)  and  A  is  self-adjoint  with  compact  resolvent.  Solutions  of  (3.1) 
are  taken  to  be  in  mild  form,  i.e., 


(3.2) 

z(t)  =  S(t)z.,  +  j  S(t  —  s)Bu(s)  ds  , 

Jo 

(3.3) 

y{t)  =  CS(t)z0  -f  C  f  S(t  —  s)Bu(s)  ds, 

Jo 

i  >  0 

We  first  consider  the  performance  index 

TOO 

(-4)  J(u)=  fMOIi;-  +  /flMf)l£H<, 

Jo 

where  y(t)  is  given  by  equation  (3.3),  and  R  >  0.  The  linear  quadratic  regulator  problem  is 
(LQR)  :  Find  u(-)  E  L2(0,oo;H)  minimizing  the  cost  functional  J  given  by  equation 


(3.4)  subject  to  the.  system  (3.2)-(3.3). 

For  the  existence  of  an  admissible  control  u  such  that  J(u)  <  oo  and  for  the  exponential 
stability  of  the  closed-loop  system  we  need  the  following  two  hypotheses. 

(HI).  The  system  (3.1)  is  stabilizable  in  the  sense  that  there  is  a  feedback  operator 
K  £  C(H,U)  such  that  the  closed  loop  semigroup  S\(t)  £  £(//)  given  by 

(3.5)  Sx(t)z  =  S(t)z  +  I  S{t  —  s)BI\SF(s)z  ds 

Jo 

for  all  t  >  0  and  z  £  H  decays  exponentially. 

(112).  The  system  (3.1)  is  detectable  in  the  sense  that  there  exists  an  operator  F  £ 
£(>',//)  such  that  the  output  injection  semigroup  SF(t)  £  C(H)  given  by 

(3.6)  SF{t)z  =  S{t)z  +  /  SF(t~  s)FCS(s)zds 

Jo 

for  all  t  >  0  and  2  £  //  decays  exponentially. 

Remark  3.1. [28,  pp.  134-135]  (i)  If  (III)  is  satisfied,  then  for  any  r0  6  //,  there  is  an 
admissible  control  uzo(-)  £  L2(0,oo;(/)  such  that  J(uzo)  <  oo. 

(ii)  Let  (H2)  be  satisfied.  Then  for  any  r0  £  H  and  u(-)  £  L2(0,oc;U)  with  J(u)  <  oo, 
z(t)  defined  by  equation  (3.2)  is  in  L2(0,  oo; //). 

Now  we  state  the  following  fundamental  results  for  (LQR)  problem,  see  [4,13]. 

Theorem  3.2.  Let  hypotheses  (HI)  and  (112)  be  satisfied.  Then  there  is  a  unique  optimal 
control  u(-)  £  L2{ 0,  oo;  U)  for  the  linear  quadratic  regulator  (LQR)  problem  and  u(-)  is  given 
by  the  feedback  law 

(3.7)  u{t)  = -R~l  B'Uz(t),  t>  0 

where  z(t)  is  the  corresponding  optimal  trajectory  and  n  £  £(//)  is  the  unique  nonnegative 
self-adjoint  solution  of  the  algebraic  Riccati  equation 

(3.8)  A'Uz  +  UAz-  \\BR~xB'l\  z  +  C‘C  z  =  0 

for  every  ^  £  T>{A).  Moreover,  the  closed-loop  semigroup  Sn(t)  £  £(//)  generated  by  the 
operator  A  —  BR~X  B"U  decays  exponentially. 
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Remark  3.3.  The  algebraic  Riccati  equation  (3.8)  can  be  defined  for  all  c  6  //,  since  the 
right  hand  side  of  equation 

(3.9)  xiaz  =  -.4*11  ~  +  nBirurwz-  cmCz 

is  well-defined  for  all  c  E  //,  T>(A)  is  densely  embedded  in  //,  and  hence  we  can  extend  the 
left  hand  side  of  equation  (3.9)  to  :  6  //  continuously. 

We  now  consider  the  control  problem,  and  its  linearization,  for  Burgers'  equation.  Define 
the  operator  At ,  as  in  Section  2,  by  Ac<f  =  c  </>"  for  all  <f  E  T>(A()  =  //2( 0,  f)  Pi  //o(0,f).  For 
the  control  input  operator  B  and  the  observation  output  operator  C  we  consider  the  Hilbert 
spaces  //  =  L2(0,  C),  U  =  R  and  Y  =  Rfc.  Assume  that  B  E  C(U.ll)  and  C  E  C(II,Y)  are 
defined  by 

(3.10)  Bu  =  b(-)u  and  Cz  =  (l(a*i  ),•••,  z{ik))< 
where  b(-)  E  //,  u  E  (/,  and  x,  E  (0,  £),  1  <  i  <  F,  are  defined  by 

(3.11)  *(*<)  =  (  *(*)  dx  . 

-bJr'-S 

In  equation  (3. 1 1 ),  S  >  0  is  chosen  so  that  (x,  —  6,  x,  +  6)  C  (0,  (!)  for  all  1  <  i.  <  k. 

Consider  the  following  lincai  control  problem. 

(LQR)a  :  Find  ii(-)  E  F2(0,  oo;  U)  minimizing  the  weighted  performance  index 


(3.12) 

roo 

J(u)=  /  {||w(0llv  +  «II«(0IIMC 

Jo 

2o(dt,  (a  >  0) 

subject  to  the  governing  equations 

(3.13) 

^z(t)  =  A'Z(t)  +  Bu(t), 

-(0)  =  -o 

(3.14) 

y(t)  =  Cz(t ),  t  >  0, 

where  /lt,  B,  C  are  as  above  and  li  >  0. 

Remark  3.4.  (i)  Equation  (3.13)  corresponds  to  the  system 
d  0^ 

—z(t,x)  =  e  -z—z(t,x)  4-  6(x)?/(/.),  t  >  0,  0  <  x  <  ( , 

at  ox 1 

r(t,0)  =  i'(t,  (!)  =  0, 

e(0,x)  =  e0(x) . 
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(ii)  For  each  i,  1  <  i  <  k,  z(ii)  given  by  equation  (3.11)  represents  an  average  value  of 
x(x)  over  a  small  neighborhood  of  x,.  We  can  regard  each  x,,  1  <  i  <  k ,  as  the  location  of 
a  “weight”. 

(iii)  The  weight  function  e2at  in  the  definition  (3.12)  of  the  cost  functional  J  will  play  an 
important  role  in  the  exponential  decay  rate,  see  Theorems  3.9  and  3.10.  However,  it  also 
gives  rise  to  the  question  of  existence  of  an  admissible  control  u(-)  such  that  J(u)  <  oo. 

For  the  control  problem  (LQR)a,  we  introduce  an  “a- shifted  ”  control  system  [15].  Let 
z(t)  =  z(t)eat ,  ii(t)  =  u(t)eat  and  y(t)  =  y(t)eat.  We  then  have  a  modified  linear  control 
problem 

(LQR)a  :  Find  u  £  L2(0,  oc;  U)  minimizing  the  cost  functional 

roo 

(3.16)  j(u)=  {||y(i)lly  +  «lkOII?r}A 

Jo 

subject  to 

(317)  ftHt)  =  (A,+aI)z(t)  +  Bu(t),  i(0)  =  .-„, 

»(«)  =  CHI),  t  >  0. 

The  solutions  for  the  system  (3.17)  again  are  taken  as  mild  solutions.  If  we  solve  the 
problem  (LQR)a  and  apply 

(3.18)  ua{t)  =  e~ai  ii(t)  (t>  0) 

to  the  original  control  system  (3.13)-(3.14),  then  the  resulting  optimal  trajectory  za(t)  will 
satisfy  the  inequality 

(3.19)  IIM0llll< 

where  M  >  1  is  a  constant  and  a  >  0  is  the  desired  degree  of  stability. 

Remark  3.5.  A  discussion  of  the  “a-shifted”  problem  for  finite  dimensional  systems  first 
appeared  in  [1].  Anderson  and  Moore  showed  that,  for  finite  dimensional  systems,  the  control 
problem  (LQR)0  is  “equivalent”  to  (LQR)a  in  the  following  senses  : 

(i)  The  minimum  value  of  J  defined  by  equation  (3.12)  is  the  same  as  the  minimum  value 
of  J  given  by  equation  (3.16). 


(ii)  If  u(t)  =  g(z(t))  is  the  optimal  control  for  (LQR);  for  some  function  g ,  then  u(t)  = 
e-atg(z(t)eat) 

is  the  optimal  control  for  (LQR)a  and  conversely. 

We  now  apply  Theorem  3.2  to  the  problem  (LQR  > in  order  to  obtain  an  optimal  control 
u  €  T2(0,  oo;  U)  for  (LQR);. 

Remark  3.6.  (i)  From  Remark  2.2,  the  spectrum  cr(Ac  -f  aJ  i  of  the  infinitesimal  generator 
A(+al  consists  of  all  eigenvalues  Aai„,  n  =  1, 2,  •  •  • ,  given  by  Aan  —  a— tn2 n2 / d2  and  for  each 
n,  n  —  1, 2,  •••,  the  eigenfunction  (f>a<n  corresponding  to  Aan  is  given  by  d>0,n(x)  =  \/2sin~x. 

(ii)  We  are  interested  in  the  stabilization  problem  for  the  system  (3.15)  with  small 
viscosity  t  =  ^  >  0,  i.e.,  high  Reynolds  number.  Let  o0  =  yr-  If  a  >  a0,  ihen  there  is  at 
least  one  positive  eigenvalue  A0il  =  o  —  75-.  Moreover,  if  a  >  o0  and  0  <  7  <  e.  then  the 
first  eigenvalue  a  -  7^-  of  A^  +  a  I  satisfies  a-:yr>a-yr>  0  and  hence  A,Toi  will 
have  at  least  one  positive  eigenvalue  for  all  0  <  7  <  e.  This  will  become  important  when 
we  apply  feedback  laws  computed  from  low  Reynolds  number  to  investigate  the  closed-loop 
response  of  the  nonlinear  Burgers’  equation  for  high  Reynolds  numbers  (see  Example  3  in 
Section  4). 

2 

Let  a  >  oq  =  7T  be  given  and  let 

2  2 

(3.20)  na  =  max  {  n  6  N  :  AQ_„  =  o  -  >  0  }. 

Since  A(  is  self-adjoint  (see  Remark  2.2)  and  the  set  :  n  =  1,2,  •••}  is  a  basis  for 

H  =  L2(0,t),  we  can  identify  z  €  H  with  the  sequence  {<  z,(f>Q<n  >}„eN.  Assume  that 
bn  E  U  and  c„  £  Y  satisfy 

OO 

(3.21)  Bu  =  {<  bn,  u  >}n€N  and  Cz  =  y~]cn  <  z,  (j>a<n  > 

n=  1 

with  mi  <  00  and  !|c„|| Y  <  00 ,  see  I28’  PP-  137-143]. 

The  following  lemma  is  an  application  of  stabilizabilitv  and  detectability  results  of 
Pritchard  and  Salamon  [28,  Section  4.2]. 

Lemma  3.7.  For  each  n  =  1,2,  ••  •  ,nQ,  let 

(3.22)  Aa.n  =  {  —  :*=1,2,  —  1}. 

n 
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Then  the  following  statements  hold. 

(a)  bn  —<  b(-),(f>ain  >7^  0  for  all  n  =  1,2,--  ,na,  if  and  only  if  the  system  (3.17)  is 
stabilizable  in  H . 

(b)  If  S  >  0  satisfies  S  <  then  for  each  n  =  1,2,- ••  ,  na,  there  exists  at  least  one 
X;,  1  <  t  <  A:,  such  that  x,  ^  Xa<n  if  and  only  if  the  system  (3.17)  is  detectable  through 

C  e  £(//,  Y). 


Proof:  (a)  From  Remark  3.6,  we  know  that  the  spectrum  cr(Ac  +  al)  of  A(  +  al  consists 
of  all  eigenvalues  Aan  =  a  —  cn2ir2 /  £2 .  Thus,  for  all  n.  n  >  na  +  1,  we  know  that  n  <  0. 
Let  Hu  be  the  linear  span  of  eigenfunctions  •  •  •  ,  <f>a,na-  Then  the  dimension  of  Hu  is  na 
and  hence  the  system  (3.17)  is  stabilizable  if  and  only  if  the  projection  of  (3.17)  onto  Hu  is 
controllable,  or  equivalently,  if  and  only  if  bn  =<  b(-),  (f)a<n  0  for  all  n,  n  =  1, 2,  ■  ■  •  ,  na. 

For  (b),  let  cn  =  (cnl,---  ,cn^),  n  =  1,2,---,  be  defined  by  equation  (3.21).  Then  we 
have 


/2£  .  rnrxi  .  nnS 


(3.23)  cni  =  — -  sin  — —  sin 

mro  £  £ 

for  1  <  i  <  k.  By  the  dual  statements  of  (a),  the  system  (3.17)  is  detectable  through 
C  €  £(//,  Y )  if  and  only  if  cn  /  0  for  all  n  —  1, 2,  •  •  •  ,  na.  Hence,  (b)  holds. 


Remark  3.8.  If  na  =  1,  then  Xa<i  is  the  empty  set. 


We  now  return  to  the  original  control  problem  (LQR)a.  The  following  theorem  is  the 
main  result  for  our  control  problem  (LQR)a. 

Theorem  3.9.  Let  a  >  ao  be  given.  Suppose  that  b(-)  £  H  =  L2( 0,£),  6  >  0,  x,,  1  <  i  < 
k ,  satisfy  the  conditions  (a)  and  (b)  of  Lemma  3.7.  Then  there  is  a  unique  optimal  control 
u0,£(-)  €  L2(0,oo;R)  for  the  problem  (LQR)a  such  that 

(3.2i)  uaAt)  =  -R~lB*na,(  za,({t),  t>  0 , 

where  cai<:(t)  is  the  corresponding  optimal  trajectory  and  no  t  €  £(//)  is  the  unique  non¬ 
negative  self-adjoint  operator  satisfying  the  algebraic  Riccati  equation 

(3.25)  {A(  +  a/)*na,f2  +  UaAA,  +  al)z  -  Ua<,BR~' B‘Uat(z  +  C'Cz  =  0 

for  every  z  £  V{At)  —  H2(0,£)  H  Hq( 0,f).  Moreover,  the  closed  loop  semigroup  *?natf(0  £ 
£(//)  satisfies  the  following  stability  property 

(3.26)  IIV,(')IU(«><Me-<“+“" 
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for  some  constants  A/  =  A/(«,  r)  >  1  and  >  0. 

Proof:  By  Lemma  3.7,  we  know  that  the  o-shifted  control  system  (3.17)  satisfies  all 
hypotheses  (111)  and  (112)  (with  c(t),  y{t),  »(/)'  A,  ^(0  and  J  replaced  by  ~(/).  y(t),  u(t ), 
At  +  al ,  S(f)  and  J,  respectively  ).  Hence,  by  Theorem  3.2,  there  is  a  unique  optimal  control 
u(t)  for  (LQR)q  and  the  corresponding  closed-loop  semigroup  S(t)  decays  exponentially,  i.e., 

(3.27)  ||5(0IU(//)<  A/c_urt  t>0 

for  some  constants  M  =  A/(a,  c)  >  1  and  ui  —  Lj(a,c)  >  0.  Moreover,  ii(t)  is  given  by 

(3.28)  ii(t)  = -R-'B'WaJ{t) 

where  z(t)  is  the  resulting  optimal  trajectory  for  the  o-shifted  system  (3.17)  and  Ll^.t  is 
the  unique  nonnegative  self-adjoint  solution  of  equation  (3.25).  Since  the  semigroup  S(t)  is 
generated  by  A(  -f  al  —  BR~1BmTJat(,  the  infinitesimal  generator  of  the  closed- loop  semigroup 
SYi0  t{t)  for  the  original  system  (3.13)  is  At  —  BR~X  BmW0^.  Hence,  5na  ,(t)  =  S(t)c~at  and, 
by  the  relation  (3.27),  Snar(t)  satisfies  the  inequality  (3.26)  with  M  =  M.  Moreover,  the 
optimal  control  u0if(t)  for  (LQR)„  is  given  by  the  formula  (3.24),  since  uQ<l(t)  =  ii(t)e~c,t  = 
— R~x B*UQ'tz(t)e~at  =  —  R~x  B*Ua^za^(l),  where  za>((t)  =  z(t)e~at  is  the  corresponding 
optimal  trajectory  for  the  original  system  (3.13).  This  completes  the  proof. 


The  optimal  control  na,{(-)  G  L2(0,  oo;R)  obtained  in  Theorem  3.9  is  given  by  the 
feedback  law  (3.24).  Note  that  uai((t)  depends  on  a  and  c.  Define  the  feedback  operator 
Ka,tZC(H,U)  by 

(3.29)  Ka,t  =  -RT'B'X\att. 

Then  the  optimal  control  f/a><(/),  t  >  0,  is  given  by 


(3.30) 


^a,<(0  Aft|(  ^or,((f) 


and  the  infinitesimal  generator  for  the  closed-loop  semigroup  Sna  {{t)  is 


(3.31)  At  +  BKa,c  =  A,  -  BR-'B-ll^. 

Recall  that  //  =  L2(0,()  and  U  =  R.  Thus,  by  Riesz  representation  theorem  (see,  e.g., 
[11,  p.  13]),  there  is  a  unique  feedback  gain  function  ka  l(-)  £  L2(0,f)  such  that 


(3.32) 


ds 
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for  all  ~  £  L2{0J). 


We  now  have  the  following  main  result  of  this  section. 

Theorem  3.10.  Let  o  >  a0  be  given.  Suppose  that  b(-)  £  II  =  L2( 0,£).  S  >  0  and 
x,,  1  <  i  <  k,  satisfy  the  conditions  (a)  and  (b)  of  Lemma  3.7.  Let  ka((-)  £  II  be  the 
linear  feedback  gain  function  defined  by  the  formula  (3.32).  Then  there  exist  constants 
p  =  p(a,t)  >  0  and  A/  =  M(a,t)  >  1  such  that  for  any  initial  data  ^o(')  G  //q(0. £),  with 
ll-oll//^  <  ^7,  the  controlled  Burgers’  equation 

—  c(Lx)  =  c  —  -  z(f, x)  —  r(f, x)  +  6(x)  /  ka^(s)z(t,s)  ds 

(3'33)  z(L0)  =  .-(<,  f)  =  0, 

2(0,x)  =  *o(*)G  //o(0,T) 

has  a  unique  (strong)  solution  and  the  solution  .:(/)(•)  =  z(t,  ■)  satisfies  the  following  stability 
property 


(3-34)  II  *(0  II//'  <2A/C-°‘||2o(-)ll//i- 

Proof:  Let  the  operators  B ,  C  and  7\a,e  be  given  by  equations  (2.4),  (3.10)  and 
(3.29).  Define  the  nonlinear  function  /  :  //q(0,  I)  — >  L2( 0,  f)  by 

(3.35)  f(z)  =  BKa,tz-zz' 

where  '  —  Then,  the  map  /  satisfies  the  hypothesis  (F)  in  Section  2,  since  for  any 
2l,  ^2  G  //<5(0,  f), 

(3.36)  1 1  /( -1  )  -  /(22)  Il/.2(0,q  <  (  ||  Rha.(\\c(H)  +  1 1^1 1 1//1  +  11-211//^  )  ||~1  -  ^2  1 1  f{i  ■ 

Note  that  the  operator  BI\ni  is  bounded  on  the  state  space  //  =  L2(0,f).  Thus,  by  Theorem 
2.1,  we  have  a  unique  local  (strong)  solution  of  equation  (3.33). 

Let  Sxnt(t) ,  t  >  0,  be  the  analytic  semigroup  on  II  generated  by  the  operator  At  + 
BKai.  Then,  by  Theorem  3.9,  S^at(t)  satisfies  the  inequality 

(3.37)  l|.V.(0*-|lH<Moe-|"+"|,|k'll» 
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for  all  z  €  H  and  for  some  constants  Ma  =  A/Q(a,e)  >  1  and  ui  =  u(a,e)  >  0.  Let 
a  <  0  <  a  -f  uj.  Then,  there  is  a  constant  Ma  =  A/0(a,e,/3)  >  1  such  that 

(3.38)  I|Sk„,<<HIh'  <M«e-0,\\z\\Hi 

(3.39)  ||  SK„,(()  2 1|„,  <  U.  -L e-«  ||.-||„ 

for  all  2  (E  //q.  Let  M  =  max{  A/„,  A/a  }  and  choose  p  >  0  with  0  <  p  <  Then  it  is 

easy  to  see  that 

r°°  1  1 

(3.40)  pM  Te-(,3-a)s  ds  < 

Jo  Vs  2 

Thus,  by  arguments  similar  to  those  in  the  proof  of  Theorem  2.7  together  with  inequalities 
(3.38)-(3.39)  and  the  expression 

(3.41)  z(t)  =  S#f0,«(0*o+  /  SKa  ((t  -  s)g(z(s))ds, 

Jo 

the  unique  global  solution  z(<)(')  =  2(L‘)  for  the  controlled  Burgers’  equation  (3.33)  exists 
and  satisfies  the  inequality  (3.34),  where  g{z(t ))  =  —  z(t)z'(t). 

4.  Approximation  and  numerical  results 

In  Section  3,  we  considered  a  linear  quadratic  regulator  problem  (LQR)^  to  obtain  a 
desired  degree  of  stability  for  the  solution  of  the  closed-loop  Burgers’  equation.  In  this 
section  we  consider  an  approximation  scheme  for  (LQR)0  and  give  some  numerical  results. 

We  first  introduce  an  abstract  approximation  scheme  for  the  problem  (LQR)  based  on 
the  results  of  Banks,  Kunish  [4]  and  Ito  [19]  and  then  apply  the  scheme  to  get  the  optimal 
control  u0it(-)  €  T2(0,oo;R).  Throughout  this  section,  we  assume  that  R  =  /,  the  identity 
operator  on  the  control  space  (7,  (_ee  Section  3).  For  more  approximation  schemes  for  the 
linear  regulator  problems,  see  [2,4,14]  and  references  given  there. 

We  formulate  a  sequence  of  approximate  regulator  problems  and  present  a  convergence 
result  for  the  corresponding  Riccati  operators.  Throughout  this  section,  we  use  superscript 
N  in  the  designation  of  subspaces,  operators  and  matrices  in  the  jV-th  approximating  system 
and  control  proolem,  like  HN ,  AN ,  BN ,  etc.  Hence  the  superscript  N  indicates  the  order  of 
approximation. 
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Let  II N ,  N  —  1,2,-  ■  ■ ,  be  a  sequence  of  finite  dimensional  linear  subspaces  of  II  and 
PiV  :  II  — »  H N  be  the  canonical  orthogonal  projections.  Assume  that  S‘s  (t)  is  a  sequence 
of  Co-semigroups  on  H N  with  infinitesimal  generators  AN  G  £(//'*).  We  then  consider  the 
family  of  regulator  problems: 

(LQR)/V  :  Minimize  Jn(zq,u)  over  u  6  L2(0,oc;f/)  subject  to  the  control  system 

(4.1)  zN{t)  =  SN(t)zg  - f  f  SN{t-s)DNu(s)ds, 

Jo 

(4.2)  yN(t)  =  CNzN(t), 

where  2^(0)  =  =  P‘N z0  and 

(4.3)  •/l'W,<*)=  /  {||/«  +  IM0llM*- 

Jo 

Remark  4.1.  If,  for  each  N,  (AN,B,W)  is  stabilizable  and  (AN ,CN)  is  detectable,  then, 
by  Theorem  3.2,  there  is  a  unique  optimal  control  uN(t)  for  the  finite  dimensional  problem 
(LQR)N  and  it  is  given  by 

(4.4)  uN(t)  =  -(BNynNs»(t)zZ, 

where  S£[(t)  is  the  Co-semigroup  on  IIN  generated  by  AN  —  BN (BN)*T[N  and  IIN  G  £(//A) 
is  the  unique  nonnegative  self-adjoint  solution  of 

(4.5)  (A^)*n'v  +  n^A"  -  n^i?;v(z?^)*n;v  +  ( cNycN  =  o. 

For  the  finite  dimensional  approximation  systems,  it  is  not  clear  that  (AN,BN)  is  sta¬ 
bilizable  even  if  the  original  system  ( A,B )  is  stabilizable.  Similarly,  it  is  not  clear  that 
the  detectability  property  of  (A,C)  is  preserved  under  the  finite  dimensional  projections. 
Another  question  we  have  to  consider  is  the  convergence  of  approximates  ITV  and  uN(t)  to 
the  infinite  dimensional  solutions  n  and  u(<),  respectively.  For  these  reasons,  we  need  the 
following  assumptions. 

Let  SN(t)  =  eAN\  t  >  0. 

(Al):  For  each  2  G  //,  SN(t)PN  z  — ►  S(t)z  and  5N (t)m PN z  — >  S(t)'z ,  where  the 
convergences  are  uniform  in  t  on  bounded  subsets  of  [0,oo). 

(A2):  (i)  For  each  n  G  U,  BNu  — *  Bu  and  for  each  2  G  //,  (BN)" PN z  — *  B‘z. 
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(ii)  For  each  c  G  //,  P‘%  z  — >  C z  and  for  each  y  G  Y,  {CN)*y  - ♦  C*y. 

(.43):  (i)  1'he  family  of  systems  (A‘W,B‘W)  is  uniformly  stabilizable,  i.e.,  there 

exists  a  sequence  of  operators  K  v  G  C(IIN,U)  such  that  sup  ||A',V||  <  oc  and 

\\e^AN+BNKN)tPs\\C{ll)  <  A/,e-u'1\  <>0, 

for  some  positive  constants  Mt  >  1  and  u.'i  >  0  which  are  independent  of  N. 

(ii)  The  family  of  pairs  (AN  ,CN)  is  uniformly  detectable,  i.e.,  there  exists 
a  sequence  of  operators  Av  G  C(Y,II N)  such  that  sup  HF^H  <  oo  and 

<  A i  >  0, 

for  some  constants  A/2  >  1  and  u?2  >  0  which  are  independent  of  Ah 

Remark  4.2.  (i)  The  condition  (2l3)(ii)  is  a  relaxation  of  the  coercivity  assumption  in  [4] 
(see  also  [19,  p . 3 ] ) . 

(ii)  Suppose  that  DN  =  PN D  and  CN  =  CPN .  Then  (A2)  holds,  since  it  follows  from 
(211)  that  PNz  ->  z  for  all  2  G  H. 

By  simple  modification  of  results  from  [19,  Theorem  2.1]  and  [4,  Theorem  2.2],  we  have 
the  following  fundamental  convergence  results. 

Theorem  4.3.  Let  (A,D)  be  stabilizable  and  {A,C)  be  detectable.  Suppose  that  (/ll)  — 
( v43 )  are  satisfied.  Then,  for  each  N,  the  finite  dimensional  algebraic  Riccati  Equation 

(4.5)  admits  a  uique  nonnegative  self-adjoint  solution  II N  such  that  sup{  | |T1',V | |,c(//jv)  :  N  = 
1 , 2,  •  ■  •  }  <00  and 

(4.6)  n *PNz  — ->  Hz 

for  every  2  G  //.  Moreover,  there  exist  positive  constants  A/3  >  1  and  u>3  (independent  of 
N)  such  that 

(4.7)  ,|e(4"-B"(*"rn <  Mzt-"3*,  t  >  0. 


For  the  uniform  stabilizability  and  detectability  assumptions  (T3)  we  introduce  a  sesquilin- 
ear  form  af(-,-)  :  V  x  V  — >  C  defined  by 


a((z,w)  —  I  tz'(x)w'(x)dx ,  z,w£V, 

Jo 
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where  V  =  //<](0,1).  Note  that,  to  allow  the  use  of  the  theory  of  sectorial  operators  and 
sesquilinear  forms  in  discussing  the  spectra  of  various  operators,  we  assume  in  defining  gc(-,  •) 
that  the  functions  in  V  are  complex  valued.  It  is  easy  to  see  that  the  sesquilinear  form  af(-,  •) 
is  I/-coercive  [36,  p.  274],  i.e., 

(4.9)  |0((2,te)|  <  e  ||~||v  |M|v  i  (continuity), 

(4.10)  Hea^Ziz)  +  ~i\\z\\2H  >  t\\z\\\, ,  (Garding’s  inequality) , 

for  all  z,w  £  V  and  7  >  e  >  0.  Furthermore,  it  follows  from  the  bounds  (4.9)  and  (4.10) 
that  there  exists,  in  a  unique  manner,  an  operator  At  £  £(K  V*)  such  that 

(4.11)  ac(z,w)  =<  —Acz,w  >v,v  and  ac(z,w)  =<  —A*w,  z  >v,v 

for  all  z,w  £  V  (see,  e.g.,  [36,  pp.  271-275]). 

Turning  next  to  specific  approximations  for  ( LQR)a ,  we  divide  the  unit  interval  [0,1]  into 
N  +  1  equal  subintervals  to  get  [x,-,  x,+1],  x,-  =  i  =  0, 1 ,  •  •  ■  ,  N.  For  each  i,  1  <  i  <  Ar, 
let  h?(x)  denote  the  linear  spline  basis  function  defined  by 

{(Af  +  l)(x  -  x,-_i),  x,-_  1  <  x  <  x, 

-(N  +  l)(x  -  x,+i),  x,  <  x  <  xi+, 

0,  otherwise. 

Let  Hn  be  the  Ar-dimensional  finite  element  space  given  by 

N 

(4.13)  Hn  =  {^af(x)  :  2,-  £  R,  i  =  1,2,  •  ■  •  ,  Af  }. 

i=i 

Then  we  have  a  sequence  of  finite  dimensional  (real)  subspaces  HN  C  V,  n  =  1,2,  •••. 
Moreover,  it  is  well-known  [31], [19,  p.  15]  that  the  family  of  HN  satisfies  the  following 
approximation  condition: 

(APP)  :  For  each  z  £  V,  there  exists  an  element  zN  £  //'v  such  that 
\\z  —  z^ 1 1 v  <  e(Ar),  where  e(N)  — *  00. 

Let  P  :  H  — +  HN  be  the  canonical  orthogonal  projection  onto  HN .  Then,  from  the 
approximation  property  (APP),  it  is  a  trivial  matter  to  see  that  PN z  — ►  z  as  A^  — +  00,  for  z  £ 
H .  For  the  finite  demensional  regulator  problem  (LQR)%  we  choose 

(4.14)  Bn  =  PnB  and  CN  =  CPN. 
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Then  the  conditions  (A2)(i),(ii)  follow  from  Remark  4.2. 

To  obtain  a  representation  A ^  of  Ae  on  consider  the  restriction  of  the  sesquilinear 
form  ae(-,  •)  to  HN  x  HN.  We  then  have  a  representation  A^  of  At  satisfying 

(4.15)  a((z,w)  =<  —A^ z,w  >  and  at(z,w)  =<  —{A1}' )*w,  z  > 

for  all  z,w  €  HN .  Equation  (4.15)  follows  from  the  fact  that  HN  is  a  real  Hilbert  space.  We 
know  also  that  A ^  =  (A^)*,  since  Ac  =  A*. 

Remark  4.4.  Since  HN  C  H,  by  equation  (4.15),  it  is  easy  to  see  that  for  any  A  £  cr(A^), 
TZe  X  <  —7  <  —  e. 

Let  SN(t)  be  the  Co-semigroup  generated  by  A^ .  Then  the  conditions  (Al)(i),(ii)  fol¬ 
low  from  the  results  of  Banks  and  Kunish  [4,  Lemma  3.2].  Note  that  SN(t)  =  (SN(t))’. 
For  the  condition  (A3)(i)  we  need  a  certain  preservation  of  exponential  stabilizability  under 
approximation  (  (POES)  in  [4]).  The  following  result  is  taken  from  [4,  Lemma  3.3]. 

Theorem  4.5.  Let  (Ac,  B)  be  (exponentially)  stabilizable.  Suppose  that  the  approximation 
condition  (APP)  holds.  Then  the  approximations  defined  through  equations  (4.14)-(4.15) 
satisfy  the  condition  (A3)(i),  i.e.,  the  family  of  pairs  (Af\  BN)  is  uniformly  stabilizable. 

By  the  dual  arguments  of  Theorem  4.5  we  can  see  that  the  condition  (A3)  (ii)  holds 
under  the  assumption  that  (Af,C)  is  detectable.  We  summarize  our  discussion  up  to  this 
point  as  the  following  theorem. 

Theorem  4.6.  Let  (Ae  -fi  a/,  B)  be  stabilizable  and  (A{  +  a/,C)  be  detectable.  Let  A 
BN ,  CN  be  defined  as  in  equations  (4.14)  and  (4.15).  Then  we  have 

(4.16)  n^PNz— .n„*,  z£H, 

and 

(4.17)  SN(t)PNz  — ♦  S(t)z,  z€  H, 

where  the  convergence  is  uniform  in  t  on  bounded  subsets  of  (0,oo),  PN  is  the  orthogonal 
projection  onto  //'v,  and  n^£  satisfies 

(4.18)  (A?  4-  o/N)Tl"  +  n ?JA?  +  aIN)  -  U^BN{BNyn^  +  (CN)*C*  =  0. 
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Remark  4.7.  Note  that  S(t)  =  S(t)e  at ,  SN(t)  =  SN(t)e  at,  where  S(t),  SN(t),  S(t)  and 
S*(t)  are  semigroups  generated  by  /1(,  ,4'v,  Ac  -f  al  and  A^  +  a7v,  respectively. 

Next,  consider  the  matrix  representations  of  operators  on  the  space  HN .  Let  the  ap¬ 
proximate  solution  zN(t,x)  of  z(t,x)  on  HN  be  given  by 

N 

(4.19)  zN(t,x)  =  Y.zAt)l>  IV) 

1=1 

for  some  z^(t)  G  R,  i  =  1,  •  •  •  ,N.  Then,  from  equations  (4.14)  and  (4.15),  we  have  a  finite 
dimensional  ODE  system 

(4.20)  [G"]|{zA'(0)  =  K'V]{2N(0)  +  {£>(0, 

where  {**(<)}  =  [**(<)»  •  ’  *  »  zn( 0]T> 

[GN]  =  [<  h?,  h?  >]/Vx^ 

4  10  0  ...  0' 

1410  ...  0 

0141  ...  0 

*  *  *  5 

0  1  4  i 

0  0  1  4  J  NxN 

'-2  1  0  0  ...  0 

1-2  1  0  ...  0 

0  1-21  ...  0 

(4.22)  [A*]  =  e(JV  +  1) 

0  1-2  1 
.0  ...  0  1-2 

(4.23)  {£"}  =  [<  6,  h?  >,  <  b,  h?  >,  •  •  •  ,  <  6,  h%  >]T, 

where  <  b,h y  >=  b(x)h^ (x)dx,  1  <  j  <  N.  Since  [G^]  is  invertible,  by  multiplying 
[GA]_1  to  both  sides  of  (4.20),  we  get 

(4.24)  |{r"(()>  =  «]{*"(()) +  {B>(1),  {;»}  = 

where 


(4.25) 


[A?]  =  [G^]-1^],  {5n}  =  [Gn]~'{Bn} 


and  {^'}  =  [Gv]  *[<  z0,  h?  >,  ■  ■  •  ,  <  z0,  A#  >]r. 


Next,  consider  a  representation  C'N  of  the  operator  C  on  IJN .  It  is  easy  to  see  that 
C'v  :  Hn  — >  R*  is  given  by 

(4.26)  [CN]  =  P(x,)],xn, 

where  Ajv(xt)  =  T  J^'-s  Ajv(x)dx,  l<i<k,\<j<N. 

Finally,  we  have  a  finite  dimensional  Riccati  equation  (4.18)  and  the  corresponding 
feedback  gain  operator  K^,  given  by 

(4.27)  A'"  =  -(B")‘ny. 

Therfore,  the  closed  loop  system  (4.24)  can  be  represented  by 

0.28)  ^U"(0}  =  (<  +  B'VA'")U'V(f)},  {*-'>))  = 


Now,  we  discuss  an  algorithm  for  finding  the  unique  nonnegative  self-adjoint  Riccati 
solution  for  equation  (4.18).  We  employ  the  Potter’s  method  [30]  to  obtain  FI^e.  The  first 
step  in  Potter’s  method  is  to  form  2 N  x  2 N  matrix 


(4.29) 


Mn 


(A?  +  aiNy  (cNycN 

BN{BNy  -(Ay  +  aIN) 


Next,  find  all  eigenvalues  and  eigenvectors  of  MN 


and  form  the  matrix 


(4.30) 


ZN  = 


N 


Qi 

q y 


where  the  columns  of  ZN  are  the  eigenvectors  of  MA  corresponding  to  the  eigenvalues 
with  positive  real  part.  When  eigenvalues  occur  in  complex  conjugate  pairs,  so  do  the 
eigenvectors.  In  this  case,  the  real  and  imaginary  part  of  the  eigenvector  each  forms  a 
column  of  Z* .  Finally,  the  solution  to  the  Riccati  equation  (4.18)  is  given  by  the  formula 

n"  =  Q?(Q?r'. 

Remark  4.8.  From  the  numerical  results  we  found  that  the  Riccati  solution  operators  Il£( 
blow  up  when  the  viscosity  e  >  0  goes  to  0  for  fixed  a  >  0.  Also,  when  o  goes  to  infinity 
with  e  fixed  the  same  phenomenon  has  been  observed. 
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Finally,  the  finite  dimensional  approximation  for  the  controlled  nonlinear  Burgers’  equa¬ 
tion  (3.33)  is  given  by 

<4.31)  =  (Wi  +  up*!,)}  +  /vu='v«))).  {=-v<»)}  = 

where  [,4^],  { }  are  defined  as  in  equation  (4.25),  and 


fS(l:Nm  =  [GN}-' m 


Vi-1  rNi  f 


7,v({=,v(0}) 


--(*iW  -  +  (-po)2 

-(-iv-2(0)2  ~  (~;v_2(0)(-;v-i(0)  +  (~.v_i(0)('\(0)  +  (c;v(/'>)2 
~(4-x(t))2  -  (^-At)Mr(t)) 


where  [G,,v]  is  defined  by  equation  (4.21).  do  solve  the  nonlinear  0 Did  system  (4.31)  we  use 
the  4-th  order  Runge-Kutta  method  [33], 

Throughout  the  rest  of  this  section,  we  discuss  how  our  results  work  for  relaxation  of 
"steep”  gradient  of  the  solution  for  Burgers’  equation  through  numerical  experiments. 


For  numerical  examples,  the  length  i  for  space  domain,  the  Reynolds  nummber.  Re,  the 
initial  function  2o(-)  G  //q (0,  1)  and  the  control  input  function  b(-)  G  T2( 0,  1)  will  be  chosen 
as  1,  60,  sin  nx  and  cx,  respectively.  Thus  the  governing  equation  is  given  by 


(4  32) 

1  1  -(TO)  =  ~(T  1 )  =  0 

r(0,  x)  =  sin  7r.r, 

where  the  feedback  gain  function  ka  t(-)  G  T2( 0, 1)  will  be  determined  by  the  desired  degree 
a  >  0  of  stability  and  the  action  of  output  operator  C. 

The  “robustness”  of  the  feedback  controller  exhibited,  for  example,  in  Figure  2.  will 
be  discussed  in  Example  3.  For  this  particular  example,  Reynolds  numbers  60,  SO,  100  and 
120  are  » hoson. 

Remark  4.9.  (i)  From  the  numerical  experiments,  we  found  that  if  Reynolds  number. 
Re.  is  less  than  60,  then  the  diffusion  phenomena  dominate  convection  phenomena.  In  this 
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case,  the  formation  of  steep  gradient  due  to  convection  term  —z(t,x)-^z(t,x)  of  the  open- 
loop  solution,  i.e.,  kat(-)  =  0  in  equation  (4.32),  is  not  clear.  But,  for  Reynolds  number 
greater  than  60,  the  open-loop  solution  creates  “sharp"  gradient  in  finite  time,  (see  Figures 
1,  10,  12  and  14).  Of  course,  the  solution  dies  out  eventually,  because  of  the  diffusion  term 
nedx>z(t’ x)- 

(ii)  The  control  input  function  b(x)  =  cx  is  defined  for  all  x  £  [0, 1].  Thus,  the 
feedback  control  acts  on  the  whole  domain  [0, 1].  But,  one  can  choose  any  L2-function 
b(-)  £  T2(0, 1)  satisfying  the  stabilizability  condition  in  Lemma  3.7.  In  fact,  b(x)  =  ex 
satisfies  the  stabilizability  condition  for  any  desired  degree  of  stability  a  >  0,  since  the 
coefficients  6„,  n  =  1,2,-  •  •  ,  representing  input  function  b(-)  are  not  zero,  i.e.,  bn  =  < 
/(■),  sin  7tx  >l2( o,i )=  Jo  ex  sin  nx  dx  0  for  all  n  =  1, 2,  •  •  • ,  (see  Lemma  3.7). 

(iii)  The  initial  function  20(x)  =  sin  nx  is  chosen  for  our  numerical  experiments.  Other 
typical  //((-functions  such  as  the  “hat  function”  defined  by 


(4.33) 


2o(*) 


2x,  x  £  [0,  f] 

— 2x  +  2,  x  £  [|,  1] 


can  be  used  for  initial  data.  But,  we  found  that  the  solution  of  Burgers’  equation  (4.32)  with 
initial  data  20(x)  replaced  by  the  hat  function  has  almost  similar  phenomena,  such  as  the 
creation  or  relaxation  of  steep  gradients,  as  those  of  sloution  with  initial  data  20(x)  =  sin  nx. 

To  show  trajectories  of  open-loop  and  closed-loop  solutions,  the  order  N  of  approximation 
is  chosen  as  N  =  32  for  both  cases.  And  the  corresponding  trajectories  from  time  t  =  0.0  to 
t  =  1.0  will  be  shown.  The  convergence  of  the  feedback  gain  functions  (•)  £  L2( 0, 1)  will 
be  shown  for  N  =  8,  16,  32,  64  and  128. 

Example  1.  (Bounded  Input/Output) 

The  observation  operator  C  £  C(L2{ 0, 1 ) ,  R3 )  for  this  example  is  given  by 


(4.34)  C(z)  =  (2(0.3),  2(0.5),  2(0.75)), 

where  5(x)  is  the  average  value  of  2(-)  £  L2( 0,1)  in  a  small  neighborhood  of  x,  x  = 
0.3,  0.5,  0.75,  and  defined  by  equation  (3.11),  z(x)  =  z(s)ds.  Here,  S  >  0  is  chosen 

so  small  that  each  open  interval  (x  —  8,  x  +  6)  is  contained  in  the  whole  domain  (0, 1).  The 
desired  degree  a  of  stability  are  chosen  0.3  and  0.6  for  Figures  3  and  5,  respectively.  For 


26 


both  cases,  na  =  max{  n  £  N  :  a  —  ^n2n2  >  0  }  =  1  and  hence  the  set  „YQii  defined  in 
Lemma  3.7  is  empty.  Thus,  all  assumptions  in  Theorem  3.9  are  satisfied. 

The  feedback  gain  functions  (•)  are  given  in  Figures  2  and  4.  From  these  plots,  it  is 
easy  to  see  that  control  action  is  concentrated  on  the  location  of  sensors.  This  phenomenon 
is  natural,  since  the  optimal  control  is  obtained  to  minimize  the  cost  functional  J  defined  by 
equation  (3.16)  whose  first  term  ||y(<)||y  =  I|C^(0IIr3  =  |^(^i ) |2,  where  x,  =  0.3,  0.5 

and  0.75  for  i  =1,2  and  3,  respectively.  The  corresponding  closed-loop  trajectories  aie 
shown  in  Figure  3  (for  a  —  0.3)  and  Figure  5  (for  a  =  0.6).  From  Figures  3  and  5,  we  can 
see  how  the  controllers  contribute  to  stabilization  of  the  steep  gradient  as  well  as  the  solution 
itself. 

Example  2.  (Identity  Output  Operator) 

For  this  example,  we  take  the  identity  operator  /  on  I2( 0, 1)  for  the  output  operator  C.  In 
this  case,  the  ouput  space  Y  is  L2{ 0, 1).  The  convergence  of  gain  functions  and  corresponding 
closed-loop  trajectories  for  a  =  0.3  and  0.6  arc  shown  in  Figures  6-7  and  8-9,  respectively. 
Since  the  observation  operator  is  identity,  this  example  gives  the  information  about  the 
maximal  control  action.  We  note  the  following  observation  concerning  the  convergence  rate 
of  gain  function.  Theoretically,  the  rate  is  O(^)  [19,  p.  15].  But,  in  this  example,  the  rate 
seems  to  be  faster  than  O(jj),  (see  Figures  2  and  4).  Another  observation  is  concerned  with 
the  location  of  maximal  control  action.  The  location  moves  to  the  left  portion  of  domain  as 
the  degree  of  stability  a  >  0  increases.  In  other  words,  we  should  put  more  action  on  the 
front  part  of  domain  to  get  a  higher  exponential  decay  rate  a  >  0.  (See  also  Figures  2  and 

4). 

Example  3.  (Robustness) 

In  this  example  we  show  the  robustness  of  the  feedback  controller  showed  in  Figure  2. 
The  feedback  controller  is  obtained  from  the  control  system  with  Re=60,  a  =  0.3,  b(x )  =  cx 
and  the  ouput  operator  C  defined  by  equation  (4.34). 

Figures  2,  10,  12  and  14  show  open-loop  trajectories  for  Reynolds  numbers  60,  80, 
100  and  120,  respectively.  The  corresponding  closed-loop  trajectories  are  shown  in  Figures 
2,  11,  13  and  15,  respectively.  The  order  N  of  approximation  is  chosen  as  iV  =  32  for 
Re=60,  80,  100  and  N  —  64  for  Re=120.  From  these  examples,  it  is  easy  to  see  that 
the  feedback  controller  obtained  for  Re=60  stabilizes  the  steep  gradient  of  the  solution  for 
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Burgers’  equation  with  various  Reynolds  numbers.  However,  we  see  that  the  sharp  gradient 
is  relaxed  slowly  as  Reynolds  number  increases,  (see  Figures  3,  11,  13  and  15).  Note  that, 
from  Theorem  3.10,  the  closed-loop  solution  z(t)  of  the  nonlinear  system  (3.33)  satisfies  the 
stability  property 

(4.35)  ||c(0||h>  <2A/e-"‘||r0(-)lk- 

Although  the  exponential  decay  rate  a  is  independent  of  Reynolds  number,  the  constant 
M  =  A/(a, Re)  depends  strongly  on  the  Reynolds  number,  Re. 

5.  Conclusion 

In  this  paper,  we  considered  a  feedback  control  problem  for  a  nonlinear  equation,  in 
particular,  Burgers’  equation.  The  method  consists  of  linearization  of  the  nonlinear  equation. 
We  used  the  linear  quadratic  regulator  (LQR)  problem  to  find  optimal  feedback  gains.  The 
linearized  equation  is  the  heat  equation.  It  was  also  proved  that,  under  appropriate  selection 
of  the  input  functions  and  “weights  location”,  the  LQR  problem  for  the  linearized  problem  is 
detectable  and  stabilizable.  We  then  analysed  a  “shifted  quadratic  cost”  to  construct  gains 
which  produce  a  fixed  decay  rate.  In  particular,  we  showed  that  the  closed-loop  system 
satisfies  the  inequality 

II  2(R0;  r0)  ||/yi  <  M(e)e~at  ||~o||//i , 

where  a  >  0  docs  not  depend  on  the  Reynolds  number,  but  A/(e)  does,  (see  Theorem  3.10). 

We  also  developed  a  numerical  scheme  for  computing  the  feedback  functional  gains. 
Several  numerical  experiments  were  performed  and  the  following  observations  were  made: 

1)  The  functional  gains  depend  strongly  on  the  “weight  location".  For  example,  if  the 
output  operator  C  is  given  by 

C(z)  =  (5(.3),5(.5),5(.75)), 

then  the  gain  function  is  concentrated  on  the  weight  location,  0.3,  0.5  and  0.75,  (Figures  2 
and  4). 

2)  The  closed-loop  nonlinear  system  is  stabilized  (as  predicted)  by  linear  feedback  laws. 
Moreover,  the  steep  gradients  (for  e  ~  0)  are  smoothed  out  by  feedback. 

3)  To  test  the  “robustness”  of  the  feedback  control  law,  one  experiment  was  performed. 

We  obtained  the  functional  gain  from  the  control  system  at  the  Reynolds  number, 


Re=60,  and  applied  it  to  the  closed-loop  system 

d  1  d  /*  i 

(5.i)  -z{t,x)  =  ——z(t,x)-z{t,x)—z(t,x)  +  cxJ  klc(s)z{t,s)ds 

at  Re=80,  100  and  120.  The  closed-loop  responses  are  shown  in  Figures  11,  13  and  15. 
Although  the  performance  was  decreased,  the  system  (5.1)  was  still  stabilized  and  smoothed 
out.  These  results  provided  some  insight  into  the  possibility  of  using  linear  feedback  laws 
for  nonlinear  distributed  parameter  systems. 
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Ft*ur*  5  Ckxwd  loop  ( fwdb*ck  by  *J#(  )) 


Figur*  6  G«io  Function  **_,(  )  (<*  0  3.  C  -  () 
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Figure  13  ('toted  L«xip  i  Re  100.  feedback  bjr  k^,(  )| 


Figure  14  Open  Loop  (M  *4  Re  1 30  *,<*)  »m*r) 


Figure  IS  doted  Loop  (>  64,  Re  1 20  feedback  b?  k^  ,(  ') 
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